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We propose a new detection method for ultrasound-modulated optical tomography, which allows 
to perform, using a CCD camera, a parallel speckle detection with an optimum shot noise sensitivity. 
Moreover,we show that making use of a spatial filter system allows us to fully filter off the speckle 
decorrelation noise. This method being confirmed by a test experiment. 



Ultrasound-modulated optical tomography is a new 
non invasive and non ionizing biological tissues imag- 
ing technique. In this technology an ultrasonic wave is 
focused into a diffusing medium that scatters the inci- 
dent optical beam. Due to the ultrasonic vibration of 
the medium, some of the diffused photons are shifted in 
frequency on an ultrasonic sideband. These are the so 
called tagged photons [l|, which can be selected to per- 
form imaging. The advantage of the method is its com- 
bination of optical contrast and ultrasonic resolution. 

_ Many groups have worked on that field. Marks et al. 
investigated the modulation of light in homogeneous 
scattering media with pulsed ultrasound. Wang et al. 
[H performed ultrasound modulated optical tomography 
in scattering media. Lev and al. made scattering me- 
dia study in the reflection configuration [J]. Wang and 
Ku || developed a frequency chirp technique to obtain 
scalable imaging resolution along the ultrasonic axis by 
a ID Fourier transform. Leveque and al. [l| performed 
parallel detection of multiple speckles on a video cam- 
era and demonstrated an improvement of the detection 
signal to noise ratio on ID images of biological tissues. 
This parallel speckle detection (PSD) is considered to be 
"so far, the most efficient technique for ultrasound mod- 
ulated optical tomography" WL, and is now extensively 
used in the field, combined aM Hi or n °t [1, EH with 
the frequency chirp technique. 

By analysing the PSD detection process, we show nev- 
ertheless that the PSD sensitivity is far from optimum. 
Moreover, as noticed by Leveque et al. PSD is sensi- 
tive to the " decorrelation of the speckle pattern, which re- 
duces the signal and increases the noise" (see also [9l[rfj|). 
In this letter, we propose to solve these two problems by 
adapting our heterodyne technique [ll| to PSD, i.e. by 
performing heterodyne PSD (or HPSD). In HPSD, the 
LO beam passes outside the sample. The LO field is 
thus much larger, and the detection sensitivity is much 
better. It is then possible to reach the optimum shot 
noise limit. On the other side, the HPSD LO beam is 
a plane wave. One can then separate the fc-space com- 
ponents of the detected field. By using a proper optical 



arrangement, which reduces the fc-space extend, we can 
fully filter off the decorrelation noise. 

Let us analyse the PSD detection process. The fo- 
cus point of the ultrasonic wave can be considered as a 
source of ultrasonic tagged photons, which are detected 
coherently by heterodyne detection. For a single pixel 
detector, the coherent detection selects the field within 
the spatial mode that fits with the pixel considered as an 
antenna. The collection efficiency can be characterized 
by its "optical etendue", defined as the product of the 
emitting area by the emission solid angle, which is the 
two-dimensional generalization of the usual Lagrange in- 
variant of geometrical optics. Due to diffraction, for one 
single mode, the etendue is about A 2 . For a N pixels 
detector, since each pixel is able to perform coherent de- 
tection within its mode, the etendue is about iVA 2 . On 
the emission side, the etendue is about nS, since each 
point of the sample external surface of area S diffuses 
photons in all outgoing directions (solid angle ~ tt). The 
collection efficiency 77, which is the ratio of the etendues 
is very low : r\ ~ NX 2 /irS. For a 2cm x 2cm x 2cm dif- 
fusing sample, rj rs 10~ 10 for 1 pixel, and 77 sa 10~ 4 for 
N = 10 6 pixels. The meaning of 77 is quite simple. Forget 
et al. 8] explained that, in order to detect the speckles 
efficiently, it is necessary to "position ... the camera ... 
to match the size of a grain (of speckle) with the size of 
a pixel" . The camera must be thus placed quite far away 
from the sample, and the photons, which are diffused by 
the sample in all directions, have a probability ~ 77 to 
reach the CCD. 

On the CCD surface, each photon is converted into a 
photo electron with a probability equal to the CCD quan- 
tum efficiency Q. In the optimal case, the heterodyne 
detection noise is related to the shot noise of the local 
oscillator. Accounting of the heterodyne gain, this noise 
corresponds to lphoto electron per pixel during the mea- 
surement time [12| (since both shot noise and heterodyne 
gain are proportional to the LO amplitude). PSD is far 
from this optimum for two reasons. Firstly, the PSD LO 
beam is obtained by amplitude modulation of the main 
laser. The noise is thus related to the total intensity (car- 
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FIG. 1: a) Experimental setup. L: main laser, MM: mov- 
ing mirror, AOM1 et AOM2: acousto optic modulators, BE: 
beam expander, M: mirror, PZT: ultrasonic transducer, S: 
sample, BS: beam splitter and CCD: CCD camera, b) Spatial 
filter system. El : beam diffused by the sample, Ea- tagged 
photon beam, Elo- LO beam, S: sample, D: diaphragm with 
rectangular aperture, O: lens and CCD: CCD camera. 



rier + sideband) , while the gain depends on the sideband 
only. Secondly, the LO beam passes through the sample 
and is diffused in many modes, which do not match with 
the mode of the detector. Neglecting absorbtion, and 
back reflection (which are also present), one looses here 
a factor r\ on the LO useful intensity. The LO intensity is 
then too small to get enough heterodyne gain for efficient 
detection. 

To improve the sensitivity, we propose to perform 
HPSD with the Fig£H a) setup. This setup is similar 
to the Toida et al. [13| one, but with a CCD camera in 
place of the mono pixel detector The main laser L 
is a A = 850nm, 20mW Newport 2010M laser followed 
by an optical isolator. The mirror MM splits the laser 
into two beams. The low intensity LO beam is shifted 
in frequency by A/ by the 2 acousto optic modulators 
AOM1 (80 MHz) and AOM2 (80 MHz +Af). It is 
expanded by BE (20 x) in order to get a plane wave (di- 
ameter 1.5 cm) larger than the CCD area. On the other 
side, the high intensity signal beam irradiates the sample 
in a 13cm width water vessel. 

The PZT ultrasonic transducer (Parametrics: f a = 
2.2MHz, diameter 35mm, focal length: 50mm) gener- 
ates an ultrasonic wave that is focused into the sample. 
The signal beam that is diffused toward the CCD (z di- 
rection) at both the optical carrier frequency (field El) 
and at the ultrasonic sidebands (field Ea) interferes with 
the LO beam (beam splitter BS) on the CCD camera. 
Accounting of the optical isolator, BS and water losses, 
the measured laser power reaching the sample is 2.5 mW. 




FIG. 2: Left hand side: 
I(k x , k y ). Right hand side 
normalized with shot noise (noise in D zone 



fc-space tagged photon intensity 
ID plot of 1{kx)- Vertical scale is 



The CCD camera (PCO Pixelfly: 1280 x 1024 pixels of 
6.7 x 6.7 iim, f c = 12.5 Hz , Ubit digital, 2.2% measured 
quantum efficiency at 850nm) records in real time the in- 
terference pattern on a PC computer. MM is adjusted to 
get an average of 2000 shots per pixels (1/2 full scale) for 
the LO beam, which remains ever much larger than the 
signal beam that is strongly attenuated and diffused by 
the sample. To measure the tagged photons field complex 
amplitude, we have chosen Af = f a + / c /4. 

According to the camera gain given by PCO (2.2 elec- 
trons for 1 LSB: low significant bit), we have measured, 
without signal beam, a noise corresponding to 1 signal 
photo electron per pixel (within 10%). Our HPSD setup 
performs thus a optimal, shot noise limited, heterodyne 
detection. 

Consider now the speckle decorrelation noise. We have 
first to notice that HPSD is less sensitive to this noise 
than PSD: the noise is the same in both cases, but the 
signal is larger in HPSD, because the heterodyne gain 
is higher. We will see now how to filter off this noise. 
The setup (Figfjji dashed rectangle) is modified as shown 
on FigJTJa. The lens O (focal length / = 250 mm) and 
the diaphragm D (~ 25 x 5 mm located in the O focal 
plane) collimate the field diffused by the sample (El and 
Ea), and reduce the fc-space extend in the k x direction. 
By this way, the speckle grains are enlarged in x and 
extended over several CCD pixels. Moreover, the LO 
beam is slightly tilted (in the x, z plane) making a angle 
9 with the z direction, so that the tagged photons versus 
LO beam interference exhibit vertical fringes (along y). 
It is then possible to separate the tagged photon signal 
(fringes) from the decorrelation noise (no fringes) by a 
simple Fourier calculation. 

To illustrate this point, we have studied the field dif- 
fused by a 3.5 cm phantom sample with a 15 Vpp (Volt 
peak to peak) 2MHz excitation of the PZT. From 4 suc- 
cessive CCD images signal (where h is propor- 
tional to I — \E\ 2 ) we have calculated the 4-phases com- 
plex signal S = (I\ — I3) + j-iji — h) where j 2 = —1. 
We have calculated S(k x ,k y ) = FFT (S(x,y)) by mak- 
ing a 1024 x 1024 truncation over the S(x, y) data mea- 
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FIG. 3: Tad /Id as a function of the sample position for 
different ultrasonic excitation: 1.4 (a), 2.5 (b) , 4.5 (c), and 
8 (d) Vpp on the PZT transductor. Laser equivalent power 
2.5mW . CCD quantum efficiency 2.2%. Measurement time 
0.96s per point. 



sured on the CCD (1280 x 1024) followed by a Fast 
Fourier Transform (FFT). The I(k x ,k y ) = \S\ 2 matrix 
(1024 x 1024) is imaged on Fig|2](left hand side) with log- 
arithmic arbitrary gray scale. The sum over the column 
I(k x ) = X)fe I((k x , k y )) is also plotted (right hand side). 
One can separate here the contributions of the product 
terms of / = E.E* (where E = E LO + E L + E A : E LO 
being the LO field, E L the field diffused at the carrier fre- 
quency and Ea the ultrasonic sideband tagged photons 
field). The tagged photon heterodyne term E l o-Ea* 
evolves fast is space (fringes). It yields to the A rect- 
angular bright zone for /, and to the A peak for I (angu- 
lar width d/f, angular offset 9). The speckle decorrela- 
tion noise corresponds to the fluctuations on the E L .E L * 
term, which evolves slowly in space (no fringes). It yields 
to the bright zone in the center of the fc-space image 
(k x w k y « 0), and to the triangular (convolution of 
2 rectangles) peak in the plot (B and C). The fluctu- 
ations of the E l o-Elo* term yield to the very narrow 
peak visible on the ID plot (k x = 0). The other terms 
give very small contributions. For example, the E L o-E L * 
and E^.Ea* terms evolve fast in time (2MHz) and are 
filter off by the CCD. As seen, a proper choice of the 9 
tilt allows us to separate in the fc-space the tagged pho- 
ton (A) and the speckle decorrelation noise (B and C) 
contributions to signal. For control purpose, we have 
considered the zone D, symmetric to A, where the shot 



noise only contributes. 

We have performed a test experiment with the diffus- 
ing PSD sample already studied in [f| (see Fig. 2 of [H). In 
the sample, a 4mm diameter vertical (x direction) cylin- 
dric black inked zone absorbing the light. The sample 
is slightly compressed in the z direction and its width is 
15mm. To get a pertinent information we have summed 
I(k x , k y ) in the A (tagged photons + shot noise) and D 
(shot noise alone) zones: 1a = J2k eA^-(^)- We have 
calculated I ad = I A — Id (tagged photons alone) , and 
plotted I ad /Id (tagged photon signal normalized with 
respect to shot noise). The FigJS^, b, c and d show the 
plots obtained by moving the sample (x = 0...20mm with 
0.5mm steps) with 1.4, 2.5, 4.5, and 8 Vpp on the PZT 
respectively. Each point corresponds to 3 successive ac- 
quisitions of 4 images (0.96s). The black inked zone near 
x = bmm is clearly seen. As seen, outside the absorbing 
zone, Tad /Id (— 0.035, 0.12, 0.37 and 1.2) is propor- 
tional to the square of the PZT applied voltage (i.e. to 
the ultrasonic power). By making A/ = / c /4, with no 
ultrasonic wave, we have measured the carrier field E L 
signal Tad /Id — 700. The maximum (8Vpp) ultrasonic 
conversion factor is then 1.7 x 10~ 3 . As seen, curves c 
and d exhibit roughly the same S/N (signal/noise) ratio. 
The technical noise is then the limiting factor. On the 
other side, on curves a and b, the tagged photon signal 
is lower and the S/N goes down. The shot noise becomes 
thus the limiting factor. With 2 zones (A and D), with 
ripi X w 2.10 5 pixels per zone, one expects for Tad /Id 
a shot noise of ±y / 2/n p i X = ±0.003, in good agreement 
with the noise observed on curve a. This shows that the 
sensitivity obtained by our technique is truly limited by 
shot noise. 

As seen above, our HPSD detection method presents 
many advantages for ultrasound-modulated tomography. 
It allows to perform parallel speckle detection of the 
ultrasound-modulated component with an optimum shot 
noise sensitivity, and to fully filter off the speckle decor- 
relation noise. At the end, many controls are possible 
on the data. One can measure , for example, both the 
ultrasound- modulated signal (zone A), the shot noise 
(zone D) and the speckle decorrelation noise (zone B 
and C). 

We thank the ESPCI group for their help (ultrasonic 
transducer, diffusing sample), and for fruitful scientific 
discussions. 
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